
function Results = print_results_theta_govs(results,dnames,rfe)

KR = length(rfe);
KD = length(dnames);

colNames = {'Estimate', 'StdError', 'pValue'};

rowNames = [cellfun(@(t)['PVEM x ',t],dnames,'UniformOutput',0),...
    cellfun(@(t)['PRI x ',t],dnames,'UniformOutput',0)];
rowNames = [{'prob_winning','incumbent_governor'},rowNames];
results([1:KR,KR+KD+(1:KR)],:) = []; % remove FEs
results = results([end,end-1,1:end-2],:); % rearrange

Results = array2table(results,'RowNames',rowNames,'VariableNames',colNames(1:size(results,2)));
Results.Properties.DimensionNames{1} = 'Coefficient';



